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Abstract. Thermal and ioniza tion structures of self-similar, magnetically-driven disk winds obtained in a compan- 
ion paper Garcia et al. (2001) are used to compute a series of jet synthetic observations. These include spatially 
resolved forbidden line emission maps, long-slit spectra, as well as line ratios. Line profiles and jet widths appear 
to be good tracers of the wind dynamics and collimation, whereas line ratios essentially trace gas excitation con- 
ditions. 

All the above diagnostics are confronted to observations of T Tauri star microjets. Convolution by the observing 
beam is shown to be essential for a meaningful test of the models. We find that jet widths and qualitative vari- 
ations in line profiles with both distance and line tracers are well reproduced. A low- velocity [Ol] component is 
also obtained, originating from the disk wind base. However, this component is too weak, predicted maximum ve- 
locities are too high and electronic and total densities are too low. Denser and slower magnetized winds, launched 
from disks with warm chromospheres, might resolve these discrepancies. 
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1. Introduction 

Jets from young protostars such as T Tauri stars remain 
an enigma. Indeed, despite tremendous efforts, both ob- 
servationally and theoretically, we do not know whether 
jets are driven by the disk alone (so called "disk winds"), 
the central protostar alone ( "stellar winds" ) or the inter- 
action zone between the disk inner edge and the proto- 
stellar magnetosphere ( "X- winds" ) . As argued in ( Garcia 



et al. 2001, hereafter Paper I) probing the ultimate en- 



ergy source powering T Tauri jels requires reliable obser- 



vational predictions from any dynamical joL model. To do 

so, the thermal and ionization states of the gas must then 
be solved along the flow. This remains to be addressed in 
a fully self-consistent way, leaving room for speculations. 
In contrast with recent approaches which assumed a 



driven disk wind heated by ambipolar diffusion. Indeed, 
of all possible heating mechanisms only ambipolar diffu- 
sion heating allows "minimal" thermal solutions, in the 
sense that the same physical process — non- vanishing cur- 
rents — is responsible for both jet dynamics and heat- 
ing. As a consequence no additional tunable parameter is 
invoked for the jet thermal description. We therefore ex- 
tended the pioneering work of paficr (1993a| ,[b|) by (1) using 
magnetically-driven (cold) jet solutions, self-consistently 
computed with the underlying accretion disk Ferreira 
(1997[); (2) a more accurate treatment of ionization us- 



constant temperature and ionization profile (Shang et al 



199S; Cabrit et al. 1999), in Paper I we solved a posteriori 



the thermal and ionization structures of a magnetically- 
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ing the Mappings Ic code and ion-neutral momentum ex- 
change rates including the thermal contribution. 

As in |5afier (1993a ), we obtained jets with a temper- 
ature plateau around 10 4 K, but ionization fractions re- 
vised downward by a factor of 10-100. This last fact is due 
to the previous omission of thermal speeds in ion-neutral 
momentum-exchange rates and to our different jet solu- 
tions. Moreover, the physical origin of the hot temperature 
plateau has been shown to represent a robust property of 
magnetically-driven disk winds heated by ambipolar dif- 
fusion. Last but not least, we showed that the calculated 
thermal and ionization states remain consistent with the 
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approximations used for jet calculations (see Appendix C, 
Paper I): (1) single fluid description of a (2) thermalized 
(all chemical species share the same temperature), (3) cold 
(negligible thermal pressure gradient) and (4) perfectly 
conducting plasma. We stress that such consistency checks 
are of an extreme physical importance and must therefore 
be made. 

In this paper, we confront the detailed model calcu- 
lations obtained in 



Paper I 



with observations of forbid- 
den line emission regions in T Tauri stars, believed to 



probe accretion-driven winds (Cohen et al. 1989; Cabrit 



et al. l [990t |Hartigan et al. 199E| , hereafter HEG95). The 



high-velocity emission has been spatially resolved in a 
dozen stars, and found to arise from c ollimated "micro- 
jets" e xtended over only a few 10 AU flSolf 1989j |Kcpncr 
et al. |1993|; |ffirth et al. 1994a|,[b] jBurrows et al. 1996|; 



Ray et al. 1996|; |Hirth et al. 1997) |Bacciotti fe Eisloffel 



1999|; [iJavalley-Fouquet et al. 200(]| ; |Dougados et al. 2000a| ; 
Bacciotti et al. 200C , and refs. therein). These microjets 
represent currently our best opportunity to test accretion- 
driven jet models: (1) At such small distances from the 
star, the intrinsic jet structure is (hopefully) less per- 
turbed by instabilities or interactions with the ambient 
medium than in the more distant jet sections accessible 
in infrared sources; (2) Direct, albeit uncertain, estimates 
for the disk accretion rate have become recently available 
from measurements of continuum excess in the stellar pho- 



tospheric spectrum (HEG9E, Gullbring et al. 1998). 

We point out that the "correct" jet model should 
pass all available observational tests, namely: microjet 
widths as measured by adaptive optics and the Hubble 
Space Telescope, long slit spectra, and integrated line pro- 
files, line fluxes and line ratios. Note that only the lat- 
ter con&traintsbased on integrated profiles were available 
to Saner (1993b ). We are now in a position to perform 
more detailed tests of the class of self-similar, magneti- 
cally driven cold disk winds, involving spatially-resolved 
diagnostics. We also stress that given the steep gradients 
in physical conditions at the base of the wind, meaningful 
comparison with observations requires proper convolution 
of maps and long-slit spectra by the observing beam. 

This article is structured as follows. Section 2 is de- 
voted to the predicted forbidden line emission diagnostics 
and their comparison with T Tauri stars observations. We 
conclude in Section 3 by summarizing the successes and 
failures of our model, and by discussing perspectives to 
resolve the discrepancies. 

2. Predicted forbidden line emission and 
comparison with T Tauri microjets 

We used three different models of self- similar magne tized 
accretion-ejection structures (MAES, Ferreira 1997 ), la- 
belled A, B and C. They are mainly characterized by a 
different ejection index 



measuring the local ejection efficiency, where M acc is the 
disk accretion rate at the cylindrical radius w. Note that 
£ = in a standard accretion disk. In our models, the ejec- 
tion index controls the jet opening, a smaller £ providing 
a larger jet radius. Model A is obtained for £ = 0.01, 
model B for £ = 0.007 and model C for £ = 0.005. A 
magnetized accretion-ejection structure is assumed to be 
settled around a protostar of M* = 0.5 M , between an 
inner and an outer edge of, respectively, zd\ — 0.07AU 
and vo c = 1AU. Outside this radius, we assume the pres- 
ence of a standard accretion disk fueling the MAES with 
mass (and possibly magnetic flux). The accretion rate is 
taken to vary between 10~*and 10 _5 M©yr _1 . We refer 



the reader to Section 2 in [Paper 1| for further details con- 
cerning the dynamical models and for a discussion of our 
choice of parameters. 

2.1. Calculation of synthetic emission maps and 
profiles 

Synthetic emission maps and long slit spectra along the 
jet axis were computed in the lines of [O i]6300, [S n]6731, 
Ha, and [Nn]6584, using the emissivities provided by 
Mappings Ic. Emission was assumed optically thin, sim- 
plifying calculations of the emergent profiles. This assump- 
tion is clearly adequate for forbidden lines, though pos- 
sibly not for Ha whose opacity is probably dominated 
by scattering of Lya photons. Our results for the latter 
line are therefore only indicative, a complete treatment of 
the line radiative transfer being outside the scope of the 
present paper. Therefore only forbidden lines will be used 
for comparison with observations of T Tauri stars, where 
the Ha line is contaminated by a strong magnetospheric 
contribution which will probably further affect the radia- 
tive transfer. 

In order to adequately sample the innermost jet re- 
gions where strong gradients in density and velocity are 
present, we used a logarithmic sampling in x = z/w and 
in wq between w\ = 0.07AU and m c — 1AU and pro- 
jected each cell onto a rectangular grid in the plane of the 
sky with 0.25 AU pixel size. Maps were convolved with 
a 28 AU (FWHM) gaussian beam corresponding to 0.2" 
at the distance of Taurus, typical of resolutions accessible 
with adaptive optics or the Hubble Space Telescope. 

To construct long slit spectra, line emission from each 
cell volume was distributed in radial velocity bins of 
3.5 km/s, and summed perpendicular to the jet axis, keep- 
ing only the projected distance from the star, d, as spa- 
tial coordinate. The result was convolved with a 70 AU 
(FWHM) gaussian along d and with a 10 km/s (FWHM) 
gaussian along the velocity axis, corresponding to typical 
resolutions achieved in ground-based long-slit spectra of 
micro jets (e.g. Solf fc Bohm 1995| ; Lavalley-Fouquct et al 
|2000| , hereafter LF2000). 



din Af acc (tn) 
din w 



(1) 



Note that qualitative fo rbidden line predic tions for 
model A were presented by |Cabrit et al. (19*99 ) using a 
highly simplified temperature and ionization prescription 



Garcia, P., et al.: Atomic T Tauri disk winds heated by ambipolar diffusion 



3 




range). It is also not overly sensitive to the detailed ioniza- 
tion and temperature structure, since we recover very sim- 



50 (00 150 

Projected distance from source (AU) 



Fig. 1. Predicted jet FWHM (gaussian fit) as a function 
of distance to the source for our models A, B, C (With 
£ = 0.01, 0.007 and 0.005) and two different beam sizes: 
28 AU and 2.8 AU, corresponding in Taurus to 0.2" and 
0.02" (VLT diffraction limit). Note the strong dependence 
on £ and the considerable bias introduced by the beam 
size. 



for the plateau. Comparison of the present accurate results 
with this preliminary work will be used to probe the sensi- 
tivity of maps and line profiles on the underlying thermal 
and ionization structure. 

2.2. Wind morphology: Images and jet widths 

Synthetic emission maps are presented in Fig. for our 3 
models and typical parameters. They successfully repro- 
duce two key features observed in microjets: (1) a strong, 
unresolved peak slightly shifted from the stellar position, 
(2) a jet-like body appearing collimated within 200 AU of 
the star. As in the X-wind model of Shu et al. (1994), this 



jet-like appearance is an "optical illusion" reflecting the 
collimation of the innermost, brightest streamlines; outer 
streamlines extend over a wider angle but have much lower 
surface brightness (cf. Fig. 1 in pabrit et al. 1999 ). Note 



that the MHD disk wind solutions considered by Safier 



(1993a) did not collimate over such scales. 

The shift of the unresolved peak from the star depends 
little on the value of £ or the chosen line, but decreases 
slightly for lower M acc , as heating then occurs sooner along 
each streamline (see Sect. 4, Paper I). It is ~ 5 AU at 
28 AU resolution fo r M acc — 10~ 6 M yr^ 1 , smaller than 
observed (15-30 AU |5olf fc Bohm 1993| ; [Hirth et al. 19944 
1997; Lavalley et al. 1997| ). However, this docs not repre 



sent a reliable test of the model as observed shifts may be 
severely biased by inaccuracy in continuum subtraction (a 
2% error in subtraction causes a ~ 4 AU error for bright 
jets; Garcia et al. 1990| ) and by possible obscuration by a 
flared disk. 

In contrast, Figures 1 and 2 show that the jet width 
provides a powerful test of the models: It depends very 
little on the chosen line tracer, inclination, jet outer ra- 
dius, and mass-accretion rate (at least within our explored 



ilar widths as in the simplified calculations of Cabrit et al 
(1999| ) , for the same wind solutions. For this class of mod- 



els, the jet FWHM depends mostly on the degree of open- 
ing of the innermost streamline, fixed by £ (Figs. ^|,[l]). Note 
that it is also severely biased by the beam size (Fig. |l]), 
an example of the importance of proper convolution for 
comparison with observations. 

Measurements of microjet widths at 0.1" - 0.2" reso- 



lution were recently presented by ( 


Burrows et al. 1996 


Ray et al. 1996; 


Dougados et al. 2000a 


). Jets with no 



clear contamination by bowshock wings give strikingly 
consistent results, namely a typical jet FWHM of 30- 
50 AU at a distance of 200 AU. Models A and B, with 
£ ~ 0.007 - 0.01, reproduce well these observations, while 
models with smaller £ < 0.007, in particular our model C, 
seem ruled out. 



2.3. Wind kinematics: Long-slit spectra and line 
profiles 

Long-slit spectra for our 3 models are presented in 
Fig. H for an inclination 60° from pole-on and M acc = 
10 _6 M Q yr~ 1 . Integrated line profiles for a wider range 
of inclinations and M acc are presented in Fig. ||. The 
predicted profiles successfully reproduce several observed 
trends in TTS: 

- Two velocity components are present: a low- velocity 
component (LVC) formed in very dense regions at the 
wind base, just where heating starts, and a high-velocity 
component (HVC), formed further up in the accelerated 
wind regions. 

- The LVC is stronger in [ O i] than in [N n] , as observed 
( |HEG95| , |Hirth et al. 1997 ): It is dominated by slowly 



rotating streamlines with low ionization fraction at vjq ~ 
1 AU. 

- The relative intensity of LVC versus HVC increases 
for lower M acc (Fig. ||), as observed (Hartigan et al. 1995): 
This occurs because the hot plateau is reached sooner fol- 
lower accretion rates, increasing emission from the wind 
base and hence the LVC. 

- The HVC is displaced further from the star than the 
LVC and the displacement is larger in [N n] than in [O i] 
(Fig. |), as observed ( [Hirth et al. 1997| ). The [Nil] profile 
peaks at the blue edge of the [Oi] li ne profile (Fig. |3 ), 
as also observed (Hartigan et al. 1995; Hirth et al. 1997). 
This difference stems from the slow increase in ionization 



fraction with distance (cf. Fig. 2 of Paper I), which favors 
Nil emission from the distant accelerated wind regions. 

- Long-slit spectra show apparent acceleration of the 
HVC, with the line profile peaking around —100 km/s 
close to the star and around —300 km/s at 150 AU. Similar 
gradients appear to be present in DG Tau and UY Aur 
QSolf fc Bohm 1993| ; [Lavalley et al. 1997| ; [Hirth et al. 1997| ) . 

The predicted profiles fail, however, in two respects: 
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Fig. 2. Predicted maps and long-slit spectra in various emission lines for Models A, B, C (£ = 0.01, 0.007 and 0.005 
left to right). Intensity maps are convolved with a 28 AU beam, and thick grey lines plot the jet transverse FWHM 
(gaussian fit) as a function of position. Long-slit spectra along the jet axis are convolved by a 70 AU x 10 km/s point 
spread function, representative of current ground-based spectro-imaging performances. Inclination is i = 60° from 
pole-on, and M acc = lO _6 M0yr _1 . A cross marks the stellar continuum position. Contours decrease by factors of 2. 



- The relative intensity of LVC versus HVC in [O i] is 
smaller than observed. The problem appears less severe 
in models with low £ < 0.005 (see Model C in Fig. ||), as 
wider streamline expansion decreases the density in the 
outer wind and the resulting HVC flux. However, such 
models predict jet widths that are too wide (see Fig. |l| and 
Sect. 2 .2 above). The weakness of the LVC in our models 



2.4. Integrated line fluxes 

Figure ^ plots the predicted integrated luminosities in the 
[Oi]A6300 and [Sn]A6731 lines for our 3 models and our 
whole range of M acc from 10~ 8 to 10~ 5 M Q yr _1 . For com- 
parison, observed line luminosities in T Tauri stars with 



the same range in accretion rates, taken from Hartigan 



reflect the low temperatures encountered over most of the 
wind base, where low-velocity gas is located. Note that in 



et al. (1995), are also plotted. 



the simplified calculations of Cabrit et al. (1999), which 
assumed T = 10 4 K above z/w ~ 0.7, the LVC was more 
prominent, illustrating the crucial influence of tempera- 
ture at the wind base on the resulting LVC intensity. 

- Predicted maximum blue radial velocities in line pro- 
files agree well with observations of the DG Tau microjet 
(e.g. [Lavalley-Fouquct et al. 2000| ; pabrit et al. 1999| ). 



However they appear too high compared with typical 
T Tauri stars, where |Vt,iuo| < 200 km/s (Hartigan et al 



1995; llirth et al. 1997), unless most stars are observed 
at inclinations > 80° from pole-on (Fig. |3|). Terminal ve- 
locities in the present class of cold disk wind models are 
~ 14ep( vcq )/V£, so even our highest £ = 0.01 (model A) 
tends to give excessive speeds. 



Two characteristics of the observations are well repro- 
duced: (1) We find a correlation between [Ol] and [Sn] 
luminosities, with a slope close to that observed in TTS 
(Fig. ^a). (2) The [Ol] line luminosity is seen to increase 
linearly with the accretion rate (Fig. p|b), a prop ortionality 
also observed in TTS (cf. Hartigan et al. 1995). 

The latter proportionality is an interesting conse- 
quence of the opposite dependence of n e and T on M acc in 
the hot plateau: As M acc increases from 10 -8 to 10 -5 M Q 
yr _1 , n e increases effectively only by a factor of 10 (re- 
call that ionization fraction decreases almost as 1/Af acc ), 
while the plateau temper ature decreases by a factor of a 
few (see Fig. 2 of Paper I ). Given the strong temperature 
dependence of [O i] collisional excitation rates qi u , the two 
effects cancel out and the emissivity (cx Q; u n e no in the 
low-density regime) becomes proportional to the oxygen 
density only, i.e. to M acc - 
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Fig. 4. Output fluxes for model A (triangles), model B 
(circles) and model C (squares). Model fluxes were inte- 
grated out to 200 AU from the star. Top: Integrated lumi- 
nosity in the [Sn]6731 line plotted against integrated lu- 
minosity in the [Oi]6300 line. Accretion rates M acc range 
from 10~ 8 to 10 _5 M Q yr _1 by factors of 10 (in the di- 
rection of increasing luminosity ). Observed T Tauri star 
forbidden line luminosities from Hartigan ct al. (1995 ) are 
shown as open circles, the arrows indicate upper limits 
on [S n]. Bottom: Integrated luminosities in the [O i]6300 
line plotted versus M acc for our models. T Tauri stars 



Fig. 3.1 Integrated line profiles behavior. Top: Effect of 



ejection index £: Model A = 0.01, solid), B (£ = 0.007, 
dashed) and C (£ = 0.005, dotted). Inclination is i = 60° 
and accretion rate M acc = lO -6 M0yr~ 1 . Middle: ef- 
fect of inclination angle for Model A and M acc = 10 -6 
Mq yr" 1 : % = 20° (dashed), i = 40° (dotted), i = 
60° (solid) and i = 80° (dash-dot). Bottom: Effect of 
accretion rate for Model A and inclination i = 60°: 
M a „. = lO^Msyr -1 (dashed), 10~ 7 M Q yr~ 1 (dotted), 

©yr" 



total [O i] luminosities are shown as circles ( Hartigan 
ct al. 1995| ); only stars with detected emission are pre- 



10 -6 M©yr _:l (solid) and 10" 5 Meyr" 1 (dash-dotted) 



Quantitatively, however, the predicted fluxes fail in 
that they are systematically too weak compared to ob- 
servations of TTS, by a factor 100 on average. We find 
L([Ol])/L Q ~ M acc /M Q yr -1 . I n this respect, we come to 
a very different conclusion from safier (1993b ), who con- 



sented. Models predict a roughly linear law L[Oi]/L0 ^ 
-Macc/M Q yr -1 , and lie typically a factor ~ 100 below ob- 
servations. 



eluded that his model was able to reproduce the range of 
forbidden line luminosities observed in TTS. There are two 
main reasons for this discrepancy: First, our lower rate co- 
efficients for H-H + momentum exchange result in a factor 
of 10 lower ionization fractions, n e values, and hence [O i] 
luminosities, compared to Safier (Paper I). Second, Safier 
considered mass-loss rates ranging up to lO _6 M yr~ 1 , 
corresponding to a mass-accretion rate ~ 10 -4 M Q yr~ 1 , 
i.e. an order of magnitude abo ve the m aximum M acc val- 
ues in TTS later deduced by HEG95 from veiling mea- 
surements. 
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In the next section, we show that the flux deficit in 
our models stems at least in part from a lower n c than 
observed in TTS, but possibly also from an insufficient 
total density in our models. 

2.5. Forbidden line ratios and total densities 

The ratios of the [Oi], [Sn], and [Nil] forbidden lines 
provide powerful diagnostics of excitation conditions in 
a stellar je t (see e.g. [Bacciotti fc Eisloffel 1999 ; Hartigan 
et al. 1994|) : [S n]6716/6731 is a decreasing function of the 



electronic density (until it reaches a "high-density limit" 
~ 0.45 when n c ^> 10 4 cm -3 , the critical density of the 
[Sn]6716 line); [Nii]/[Oi] is an increasing function of the 
electronic fraction / G ; and [Sii]/[Oi] decreases with in- 
creasing temperature^]. Hence, these line ratios provide 
a test of the heating/ionization mechanism in the flow, 
though they do not test the underlying dynamical solu- 
tion. 



In particular, LF2000 and Lavallcy (2000) showed how 
diagnostic diagrams based on these line ratios (see their 
Fig. 3) could be used to separate various proposed heating 
mechanisms, such as shocks, mixing-layers, and ambipo- 
lar diffusion heated MHD disk winds. For the latter case, 
LF2000| used preliminary results from the present work. 
We plot in Fig. || our complete set of results using the 
same line ratio diagrams. 

Each curve in Fig. corresponding to a given wind 
model and mass accretion rate, shows the evolution of line 
ratios with distance from the star from AU to 200 AU, 
after convolution of jet images by a 70 AU beam and sum- 



mation perpendicular to the jet axis (note that LF2000 



plotted ratios out to a distance of 600 AU, for comparison 
with observations of DG Tau, and used calculations with 
no depletion on dust, which affects the [Sii]/[Oi] ratio). 

Spatially resolved line ratios within 200 AU of the 
star are available for only four TTS microjets so far: 



HH 30 Th 28, DG Tau, and RW Au r flBacciotti fc Eisloffel 
1999|; lLavalley-Fouquet et al. 2000|; |Bacciotti et al. 2000| ; 



Bacciotti 2000| ; Dougados et al. 20001 f ). Comparison with 
our predictions, at the same beam resolution and distance 
from the star, indicates that our models have slightly too 
high temperatures but insufficient beam-averaged elec- 
tronic density and fractional ionization. In particular, the 
observed [S n]6716/6731 towards the star is always in the 
high-density limit ~ 0.45, indicating n c >2x 10 4 cm~ 3 , 
while our models predict higher ratios at the stellar po- 
sition corresponding to a beam-averaged n c ~ 0.6 — 2 x 
10 3 cm' 3 , This difference of at least a factor 10-30 read- 
ily explains a substantial fraction of the flux deficit noted 
above in Section 2.4. Two further important remarks are 
in order: 

(1) We stress that, because of volume effects and strong 
density gradients, beam-averaged line ratios differ greatly 
from local line ratios. In particular they give very little 



The latter two ratios also decrease with increasing n e when 
n e exceeds the critical density of the line in the numerator 
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Fig. 5. Line ratio diagrams for Model A (solid), B 
(dashed) and C (dotted). Each curve refers to a given 
M acc , and connects line ratios at various distances along 
the jet, from z = —50 AU to 200 AU (after convolution 
by a 70 AU beam and summation perpendicular to the 
jet axis). Distance increases with [S n]6716/6731 and with 
[N n]/[0 1] (bottom to top in both figures). Accretion rates 
vary from M acc — 10~ 8 to 10~ 5 M^yr -1 in factors of 10 
from left to right in both figures. 



weight to the inner densest parts of the wind, where the 
[S n]6716/6731 ratio is locally in the high-density limit. 
Hence it is essential to apply proper convolution by the 
observing beam to perform meaningful comparison with 
observed ratios. 

(2) All predicted line ratios, in particular the highly 
ionization sensitive [Nii]/[Oi], would agree much better 
with observations if one includes extra mechanical en- 
ergy deposition in shocks 



as demonstrated in LF2000 
(see also Dougados et al. 20001 ; Lavalley 200C| ) , although 
this would involve additional free parameters (e.g. ratio 
of shock speed to wind speed, magnetic field limiting the 
postshock compression) . 

For the latter reason, as noted by [Bacciotti fc Eisloffcj 
(1999D, the total density nn = n c / f c is a more fundamen- 



tal quantity that can be directly compared with MHD 
model predictions independently of the ionization process 
(although 7Jh can still be affected by shock compression 
and by beam averaging). Out of the 4 above-mentioned 
T Tauri stars with spatially resolved line ratios, only two 
(DG Tau and RW Aur) have esti mates of disk accretion 
rates, both of order 10~ 6 Moyr" 1 ( |HEG95|) . Observed val- 



ues of 71h as a function of distance (LF2000; ppugados 
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et al. 2000b) are 10-100 times higher than total densities 
on the innermost streamline for our models with the same 
A/acc (see Fig. 1 of Paper I). While shock compression 
could account for a factor of 10 Hartigan et al. (1994), it 
seems that our models may still have too low density close 
to the star. 

As a word of caution, we note that the above two stars 
are the strongest forbidden line emitters in the HEG95 
sample and may therefore not be representative of wind 



surfaces seem to be ruled out by observations of T Tauri 
microjets. In the context of magnetized accretion-ejection 
models considered here, a possible improvement would be 
to relax the assumption of isothermal magnetic surfaces. 
Indeed, very efficient ejection (namely £ up to 0.5) can 
be attained by allowing for some entropy injection at the 
disk surface Gasse fc Ferreira (2000| ). Such models are very 
promising for they naturally provide both higher jet den- 
sities and lower terminal velocities (see Figure 8 in Cassc 



physical conditions in an average TTS. More sensitive ob- fe Ferreira (2000[)). In this respect, we independently sup 



servations would be needed to probe n c and total density 
riH in objects of weaker wind emission and lower accretion 
rates. 



port Kwan (1997| ) 's point of view that active disk chromo- 
spheres may be necessary ingredients in T Tauri circum- 
stellar disks. 



3. Summary and Conclusion 



In a companion paper Paper I, we computed the ther- 
mal and ionization structures of self-similar, magnetically- 
driven disk winds heated by ambipolar diffusion. In this 
paper, we confronted our model predictions in forbidden 
lines with all existing observational constraints in T Tauri 
microjets, namely: (1) jet images; (2) jet widths; (3) long- 
slit spectra and integrated line profiles; (4) line fluxes; and 
(5) line ratios and total densities. Jet widths, line pro- 
files, and total densities provide powerful tests of the jet 
dynamics and collimation, whereas line ratios and fluxes 
are strongly influenced by the underlying heating assump- 
tions. Given the steep gradients in physical conditions in 
the jet, beam convolution is found to have a crucial effect 
on model predictions, and we stress that it is an essential 
step for a meaningful comparison with observations. 

The model reproduces several observed trends: (1) 
Images show an unresolved peak and an extended high- 
velocity jet, of width compatible with observations; (2) 
Line profiles present a low- velocity component (LVC), 
compact and near the star, and a high-velocity compo- 
nent (HVC), tracing the jet. These components are not 
an inclination effect but intrinsic to the jet dynamics. (3) 
The systematic differences between [O i] and [N n] profiles 
are reproduced; (4) Line fluxes are proportional to M acc ; 
and (5) the [O i] - [S n] correlation slope is recovered. 

The model also encounters several problems: (1) The 
line ratios ([S n]6716,6731, [Nii]/[Oi]) are not well repro- 
duced, and indicate insufficient ionization in the model; 
(2) Line fluxes are too low by a factor of 100, and electronic 
densities by a factor of at least 10-30; (3) The relative in- 
tensity of the low- velocity component (LVC), present in 
our profiles and formed in the inner wind regions just 
above the disk, is insufficient; (4) Wind terminal veloc- 
ities appear too high, and total densities at the wind base 
appear too small, when compared with current estimates 
in 2 bright microjets. 

Both line intensity and line ratio problems can be read- 
ily solved b y introducing an extra ionizin g mechanism, e.g . 
in sho cks ( Lavalley-Fouquct et al. 2000 ; Dougados et al 
2000b ). However, terminal wind velocities and total den- 
sities would not be affected. As a consequence, cold self- 
similar MHD disk wind models with isothermal magnetic 
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